Quantum fluctuation driven first order phase transition in weak ferromagnetic metals 
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In a local Fermi liquid (LFL), we show that there is a line of weak first order phase transitions 
between the ferromagnetic and paramagnetic phases due to purely quantum fluctuations. We predict 
that an instability towards superconductivity is only possible in the ferromagnetic state. At T = 
we find a point on the phase diagram where all three phases meet and we call this a quantum 
triple point (QTP). A simple application of the Gibbs phase rule shows that only these three phases 
can meet at the QTP. This provides a natural explanation of the absence of superconductivity at 
this point coming from the paramagnetic side of the phase diagram, as observed in the recently 
discovered ferromagnetic superconductor, UGe^. 
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The study of weak itinerant ferromagnetism, both ex- 
perimentally and theoretically, is an extremely important 
topic in understanding strongly correlated electron sys- 
tems. It has been shown over the past five years or so 
that materials which can be considered weak ferromag- 
nets (small Curie temperature) display a very wide as- 
sortment of complex phenomena and novel physical prop- 
erties. For example, UGe2, ZrZri2, and URhGe have 
been observed to be superconducting and ferromagnetic, 
whereas it had previously been expected (yet not ob- 
served) to be superconducting on only the paramagnetic 
side of the phase transition lil |2j. yf It is even more in- 
teresting because this superconductivity might be BCS- 
like, i.e., singlet pairing, as opposed to all of the recent 
models which predict triplet pairing. The s-wave singlet 
model has been considered by these authors and others 
and there are no definite answers as of yet 0, IMS 00- A 
new result reported recently yj shows that the magnetic 
transitions in the heavy fermion itinerant ferromagnetic 
superconductor XJGe-x are of first order, and therefore 
there does not exist a quantum critical point as previ- 
ously thought. 

In this Letter we propose an explanation of the ob- 
served first order magnetic transition and superconduct- 
ing behavior based on the 'Induced Interaction Model' 
first proposed by Babu and Brown [lfj and the general 
properties of a Local Fermi Liquid. As explained below, 
this analysis leads to a thermodynamically consistent, 
first order phase transition from the ferromagnetic state 
to the paramagnetic state. We end by considering some 
aspects of the superconductivity. 

The local Fermi liquid (LFL) was a concept proposed 
by Engelbrecht and Bedell to look at normal param- 
agnetic metals ^lj- It i s a generalization of the LFL 
proposed by Nozieres while studying the single impurity 
Kondo problem [l^. Blagoev et al recently studied a 
LFL to explain weak ferromagnetic metals @, . This 
was shown to reproduce non-trivial results and even the 
possibility of superconductivity. The superconductivity 
was predicted to be s-wave on the ferromagnetic side, due 



to the constraint of the LFL, and it had been speculated 
to be p-wave on the paramagnetic side. While the nature 
of the superconducting order parameter is still in ques- 
tion, the real mystery is why superconductivity is only 
found in the ferromagnetic state. 

To gain insight into this problem, we start with a sim- 
ple model that can describe a strongly correlated Fermi 
liquid, the LFL of ref.0|. This theory makes the as- 
sumption that the quasiparticle self-energy S(w) is mo- 
mentum independent. This leads to a further simplifi- 
cation in the theory since only the s-wave Fermi liquid 
parameters and scattering amplitudes are nonzero. In the 
limit of small magnetic moment, the scattering amplitude 
in the LFL can be expressed as A™ = Aq + A^a ■ a', 
where the A's are the scattering amplitudes, related to 
the Landau parameters by Aq a = Fq 0, /(l + Fq' 11 ). Note 
that here and elsewhere, the capital letter quantities have 
been made dimensionless by multiplication by the den- 
sity of states, A = N(0)a and F = N(0)f. In the 
LFL the forward scattering sum rule, which is a conse- 
quence of the Pauli principle, imposes the constraint that 
A^ = Ag+Ag = 0. It can be shown that when these con- 
straints are applied to a paramagnetic Fermi liquid, the 
system is stable against a transition to a ferromagnetic 
state and also against phase separation i.e., as F S 
gets large then F a saturates to -| and superconductivity 
in both the s and p-wave channels is suppressed. 

These results change dramatically when the LFL is ap- 
plied to a weak ferromagnetic system (see refs.0,0 for 
more details). In the vicinity of the phase transition, 
Fo a — > — 1" and Fq s — > — 1 + and both scattering ampli- 
tudes diverge, indicating an instability in the spin (Aq 11 ) 
and in the charge (^o s ) sector. The ferromagnetism re- 
mains while the singlet scattering amplitude, Aq 1 " 9 , is 
attractive. This opens up the possibility of the existence 
of s-wave superconductivity since the triplet scattering 
amplitude is strictly zero in the LFL. 

For small magnetization and low temperatures and en- 
ergies, one can see that this description remains valid 
within weak ferromagnetic Fermi liquid theory. However, 



FIG. 1: The schematic integral equations for the Landau /- 
function and the scattering amplitudes a. The fPs are the 
incoming and outgoing particle momenta. Note how the in- 
teraction in a) is graphically shown to be decomposed into 
a direct term d and the induced term. Part b) is the fully 
reducible set of diagrams. 

as mo (magnetization) gets even smaller as one moves to- 
ward the critical transition, the situation changes. Since 
in the local limit, m*/m ~ z" 1 ^ logmo, around the 
critical point the effective mass diverges, the quasiparti- 
cle residue z goes to zero, and the validity of Fermi liquid 
theory becomes questionable. Hence, the ferromagnetic 
and paramagnetic LFL states are not continuously con- 
nected through the critical point within this theory. But, 
and this is a crucial point, if the transition is 'preempted' 
by a first-order transition as to restrict or limit the diver- 
gences, then we can recover a consistent theory. This is 
indeed what is seen to happen as will be explained below. 

We have shown that the LFL yields unexpected and 
quite distinct predictions about the behavior of the para- 
magnetic and ferromagnetic Fermi liquids. Moreover, 
we argued that they are not connected by a continu- 
ous second order phase transition. What we will show 
here is that they are in fact connected by a first or- 
der transition using the 'Induced Interaction Model' 01 > 
which was further developed by Bedell and collabora- 
tors 0, This a model for self-consistently calcu- 
lating the quasiparticle scattering amplitude (fully re- 
ducible interaction) in terms of the three interaction 
channels: particlc-holc, exchange particle-hole(induced 
interaction), and the particle-particle channel. 

Including all three channels gives rise to a properly an- 
tisymmetrized scattering amplitude, a, where N(0)a = 
A. The diagrammatic structure of these equations is 
shown in Fig^ The direct interaction, d, is an anti- 
symmetrized effective two-body potential in the particle- 
particle channel. It is chosen specifically for a certain 
physical model, i.e., it contains information about the 
underlying Hamiltonian. The second set of diagrams on 
the right side of Fig^) is the exchange of topologically 
equivalent diagrams in Fig^j). Thus, the induced in- 
teraction is a purely quantum effect, arising from the 



FIG. 2: The two sets of solutions for the paramagnetic state 
(top two lines) and the ferromagnetic state (bottom two lines). 
The limiting values for large N(0)U of Fq' s are given in the 
text. Please note the change in scale of the y-axis shown for 
clarity. 

exchange diagrams that are required to antisymmctrize 
the effective two-body scattering amplitude. 

To determine the Fermi liquid parameters, we study 
the equations of Fig^ in the limit of the momentum 
transfer, q = pi — p§ = 0, In the local limit 

of these equations, the induced interactions are equiv- 
alent to the limit of the exchange momentum transfer 
(f = px — p~4 = 0. The full momentum dependence of 
the interactions on the Fermi surface has been inves- 
tigated extensively in the paramagnetic phase 0, Hif . 
and we are currently extending this approach to the fer- 
romagnetic phase |l5j . Including the full momentum de- 
pendence will not change the results we describe here in 
any qualitative way. Thus for our purposes, we will fo- 
cus here on the local limit of the model to calculate the 
quasiparticle interactions for both the ferromagnetic and 
paramagnetic LFL. In the local limit, the coupled inte- 
gral equations [HQ ofFigHt) andD}}) reduce to two 
coupled algebraic equations, 

*o = + \W + \FSA a (1) 
FS = D a + If%AZ - \FgA%. (2) 

The LFL picture, now coupled with the induced inter- 
action model, gives a description of a weak ferromagnetic 
Fermi liquid and its first order transition to the param- 
agnetic Fermi liquid. We will take the antisymmetrized 
direct interaction to be = -Dft = N(0)U/2, where U 
is the 'on-site' contact interaction, such as in the Hub- 
bard model. This is our model-dependent parameter. 

We solve eqtn's (11121) self-consistently using the above 
form of the direct interaction. We show the results graph- 
ically in FigI3 The different branches for each solution 
are described in the caption. The important consequence 
is that in the large N(0)U limit, the solutions remark- 
ably yield exactly the same results as those of the LFL, 
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FIG. 3: The shape of the chemical potential curves expanded 
around the critical U c - Note how the higher U regime is the 
one in which the ferromagnetic state has the lower energy. 
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Investigating the model more closely we will now em- 
ploy certain aspects of a spin-polarized Fermi liquid to 
the thermodynamics of the system. In this theory the 
Landau parameters are modified due to a finite magne- 
tization in the ferromagnetic state. For full details see 
[HI Kl\ and references therein. This has consequences 
for the effective masses to* of each spin species, as well 
as for the Fermi momenta k" F of each Fermi surface, since 
the magnetization dependence of these quantities is ex- 
tremely important. 

Since the multiple solutions of this model are found 
for every U at large enough U , one can imagine the sys- 
tem, for some specific U , jumping from one solution to 
the other (see Fig|2Jl. This can be shown by examining 
the chemical potential as a function of U. We can ex- 
pand the chemical potential around a certain U c which 
determines where one state gives way to another state 
of lower energy. The point at which the chemical po- 
tentials cross is the point of the first order phase transi- 
tion. To see this, we calculate the change in the chemi- 
cal potential due to the change in the magnetization for 
fixed density, n, which is given by Sfi — j(C'' — C^)6m, 
where C Ga = l/N a (0) + f aa and A CT (0) = k F m* a /2n 2 is 
the density of states at the Fermi surface of spin a (the 
tilde distinguishes the Landau parameters in the polar- 
ized state). Then the chemical potentials for the two 
phases near U c are written as follows: 

Mu) « Mu c ) + (u- u c ) dMu d = Uc) (3) 

MU) « Mc(t/c) + (U- U c ) dMU d = Uc) . (4) 

We differentiate the chemical potential in the ferromag- 
netic state implicitly through the magnetization which is 
itself a function of U. This is seen, for example, through 
the relation for the equilibrium magnetization in a weak 
ferromagnet: too ~ |1 + F§ (U)\ a , where a depends on 
the order in which the Ginzburg-Landau type expansion 



FIG. 4: The equilibrium magnetization plotted as a function 
of the effective interaction U = N(0)U for different T. The 
temperatures in this model are scaled by the spin fluctuation 
temperature, which is about 1/100 of the Fermi temperature 
for the interaction strengths of interest. The moment drops 
discontinuously to zero at the critical U c = 150 at zero tem- 
perature and for T < T c (shown by the vertical line), yet 
continuously to zero at the critical temperature T c = 0.2, 
indicating a crossover from first order to second order. See 
FigEl 



is carried to in the magnetization (a — 1/2 for maximal 
to term, 1/4 for to 6 , etc.). This holds for small Too |5|. 
It is seen in Fig|3| that for smaller U the paramagnetic 
state is favored but at the critical value U c the chemical 
potentials cross and the ferromagnetic state is the lower 
energy state. 

Extending these calculations to low but finite tempera- 
tures can be done using certain thermodynamic Maxwell 
relations. Doing this we can map out the temperature 
phase diagram. The concern will be with the chemi- 
cal potential and the pressure at finite temperatures and 
small magnetizations where we apply the same analy- 
sis we did at zero temperature. The details will be 
shown elsewhere, but the first step is to integrate the 

Maxwell relation, — {§^) Tm — (fr) > w ^ res P ec t to 

T, where the entropy density s is given by the usual low 
temperature Fermi liquid approximation s(n,m,T) = 
7(n, m)T |lfj| . This results in a magnetization-dependent 
chemical potential expansion in to and up to second or- 
der in T . A similar method is used to develop a free 
energy expansion in the magnetization and the tempera- 
ture. By differentiation of the free energy, an expression 
for the temperature-dependent magnetization can be de- 
rived. And finally the pressure, P = —f + fin + Hm, can 
be calculated from the free energy to give an expression 
in terms of small magnetization and low temperatures: 

" 2 dN n 2 ^ 2 



P(m,T) = P(0,0)+No Y T<-n ^ ( . 



T 



2 2 
+ GlTO 2 +G 2 TO 4 + G 3 TO 2 ^-T 2 + G4TO 4 yr 2 . (5) 

The coefficients Gi depend on the Landau interaction 
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FIG. 5: The T c versus P schematic phase diagram generated 
from this model. The double line indicates the line of first 
order transitions, which ends at a finite T, estimated in the 
text below. The SC dome, calculated from LFL theory only, 
is taken from ref 's [5J, |6| . The superconducting transition at 
low temperatures near the QTP is still under investigation as 
to its order. For further discussion on this point, see [Tsj . 

functions and the polarization expansion coefficients of 
quantities such as the effective mass [l5|. 

Now we can determine what happens at temperatures 
away from zero. At zero temperature the magnetization 
jumps discontinuously to zero at a certain U indicating a 
first order transition. However, as shown in FigQ] when 
the temperature is turned on, we see that for a certain 
critical U c , and above a certain temperature T c ^ 0, the 
magnetization goes continuously to zero, indicating a sec- 
ond order transition. This implies that there is a line of 
first order transitions that ends and soon becomes a line 
of second order transitions at finite temperatures. 

To 'translate' from U back to a physical control param- 
eter, such as pressure, we take advantage of the thermo- 
dynamics described above to develop the phase diagram. 
The results are shown in FigEl The steep slope of the 
transition near the critical pressure is due to the fact that 
the latent heat of the transition, which is zero at T = 0, 
is also zero or nearly zero at small temperatures. This is 
because the entropy difference between the FM and PM 
states is proportional to the DOS difference on each side, 
which to leading order is zero, and to the next higher or- 
der is 0(m 2 ), where m is small. Thus these transitions 
are truly weakly first order. What we see strikingly re- 
sembles the experimental phase diagram of UGe^- Look- 
ing at the physical values of this system, our effective cho- 
sen U c = 150. An approximate calculation of the DOS 
of UGe-i at the Fermi surface yields a value of w 20/eV 
[l9| . This puts U in a range of lOel^, a typical value for 
correlated electron systems. Our scaled T c = 0.2 gives, 
when a typical Fermi temperature of 1 — lOeV is used, a 
critical temperature of around lmeV, or about 10K, the 
same order of magnitude as seen in the crossover regime 
of UGe 2 - 

One final thermo dynamical observation can be made 
when looking at what we call the quantum triple point, 



QTP. At this point, as is shown in Fig[5j three phases 
end at zero temperature. According to the Gibbs phase 
rule [lOl ! a single component system, which we have here, 
can only accomodate a maximum of three phases coex- 
isting at a point. This restrictive condition explains why 
superconductivity can only be observed on one side of 
the QTP, and, as shown in previous studies 0, El , this 
must be the ferromagnetic side, where consequently only 
pairing in the singlet channel is attractive. This is seen 
experimentally. 

In summary, we explicitly have here a microscopic 
model that unmistakably yields a first order phase tran- 
sition from the ferromagnetic to the paramagnetic state 
by the inclusion of the quantum fluctuations that arise 
from the induced interactions. This can be considered 
the quantum analogue to the case of first order tran- 
sitions driven by classical fluctuations in certain liquid 
crystals studied by Brazovskii [2]]]. Consequently, if we 
turn off the (quantum) induced terms in our picture, the 
model reduces to a Stoner model which is just a standard 
second order transition between the ferromagnetic and 
paramagnetic phases. We show that the superconductiv- 
ity is s-wave and only exists in the ferromagnetic state, 
and that, as a consequence of the Gibbs phase rule, there 
are only three phases that meet precisely at the quantum 
triple point. 
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